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We compute the nonleptonic weak decay width of a heavy-light meson in 1+1 
spacetime dimensions with a large number of QCD colors (the 't Hooft model) 
as a function of the heavy quark mass. In this limit, QCD is exactly soluble, 
and decay modes are dominated by two-particle final states. We compare 
the results to the tree-level partonic decay width of the heavy quark in order 
to test quark-hadron duality in this universe. We find that this duality is 
surprisingly well satisfied in the heavy quark limit, in that the difference 
between the sum of exclusive partial widths and the tree-level partonic width 
approaches a constant as M ^ oo, and the deviation is well- fit by a small 
1/M correction. We comment on the meaning of this conclusion and its 
implications for the use of quark-hadron duality in hadronic physics. 
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I. INTRODUCTION 



Quark-hadron duality, in its most general form, is the notion that certain rates for 
processes involving hadrons can be computed simply as the underlying partonic rates [|l|. 
Duality allows us to compute many quantities which would otherwise be hopelessly difficult. 
One common application of duality is to the nonleptonic weak decays of heavy hadrons. 
The lore is that, for large enough heavy quark mass, duality holds in the computation of 
the hadronic width. 

Several discrepancies between theory and experiment that have recently received at- 
tention rely on quark-hadron duality. Among them are the significant difference between 
lifetimes of beauty baryons and mesons , the overestimates of the i?-meson semileptonic 
branching fraction ^ and the average number of charm quarks per B decay [Q. Because 
the limit of experimental knowledge about nonleptonic B decays is rapidly expanding, such 
issues are of great topical interest. 

But when is duality valid? In many cases duality follows from the Operator Product 
Expansion (OPE). This is the case, for example, for the rate of e^e~ —>■ hadrons and for the 
semileptonic decay rates of heavy hadrons. However, duality is applied in many other cases, 
such as in hadronic widths of heavy hadrons, for which there is no OPE. 

Reference proposes an OPE-like expansion in inverse powers of the heavy quark 
mass M, which not only incorporates quark-hadron duality as the lowest term in the ex- 
pansion, but also organizes the corrections by inverse powers of M. A main result of that 
work is the claim that corrections first appear at order The question above can be 

reformulated as, "Is an OPE-like expansion like that of Ref. valid?" 

To investigate the validity of duality it is convenient to work with a soluble model of 
strong interactions formulated as a full-fiedged field theory, so that one may test duality both 
in cases with and without an OPE. The 't Hooft model [Q, large- A^^c QCD in 1 + 1 dimensions, 
is a good laboratory for this purpose. It contains an infinite spectrum of mesons composed of 
confined quarks, realizes asymptotic freedom trivially, and inherits all the phenomenological 
consequences of large- A^^c QCD common to our universe, such as dominance of scattering 
amplitudes with the minimum number of meson states, OZI suppression, the absence of 
exotics, and others. For processes with an OPE, duality in the 't Hooft model has been 
checked explicitly However, little is known about duality for non-OPE processes. 

The reason is that, in precisely those cases for which an OPE is lacking, there is no simple 
analytical method of verifying duality, and one must resort to arguments based on numerical 
solutions. 

In this paper we compute the hadronic weak decay width r(M) of heavy "5" mesons in 
the 't Hooft model as a function of the heavy "6" quark mass M; the meaning of "heavy" 
is made precise in Sec. We compare this to the partonic (perturbative) decay width of 
the heavy quark, TpartiM), which we compute analytically. For large M we find that both 
rpart(^) and r(M) are essentially linear in M. The difference between the two appears to 
be asymptotically constant and small, indicating a small 1/M correction to the naive duality 
limit. As M increases, new hadronic decay channels become accessible, and at each of these 
thresholds there is a singular peak in r(M). Averaging r(M) over a region in M that 
includes many resonances removes these peaks but does not change the leading dependence 
on M. Our conclusion is that duality holds to leading order in M, hut unlike the OPE-like 
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expansion of Ref. /Q/, appears to have 1/M corrections. 

In Ref. |]lO|| it is argued that there is strong experimental evidence for the failure of 



duality. What is meant there is that the pattern of corrections in powers of 1/M of Ref. 
is not supported by experiment. This agrees with our result, which indicates a violation to 
duality at first order in the 1/M expansion of Ref. [^] . 

While we cannot carry over our quantitative results to the physical world of non-planar 
QCD in 3 + 1 dimensions, we believe that there is nothing intrinsic to 1+1 dimensions that 
would make duality work differently than in 3+1. The operator analysis that leads to the 
1/M expansion proceeds in 1+1 much as in 3+1. 

The paper is organized as follows. In Sec. 0, we briefly review the 't Hooft model 
and a standard method for its numerical solution. Section |T| compares features of 1+1 
dimensions, such as the nature of phase space and spin, to those in 3+1 dimensions. In 
Sec. we present the algebraic results of the inclusive parton-level calculation of widths. 
In Sec. |V|, we present the results of the exclusive calculation in the 't Hooft model. Section |V| 
gives our numerical results and a discussion of their implications, and Sec. [VI1| concludes. 

There are other nonperturbative questions of phenomenological interest in the area of 
hadronic B weak decays for which there is an established lore. One can test any of these 
hypotheses in the 't Hooft model. Of particular interest is the notion that contributions 
to decay amplitudes from different underlying quark diagram topologies contribute with 
distinct weights, which can very much suppress the amplitude from a given topology. For 
example, the "annihilation" diagrams, in which the valence quark-antiquark pair annihilate 
through a weak current, are supposedly suppressed relative the "spectator" diagrams by a 
factor of Jb/Mb, where Jb is the B decay constant and Mb its mass. We will address this 



question in a separate publication [11 



II. THE 'T HOOFT MODEL 

The success of 't Hooft's method of solving a strongly-coupled theory rests on two as- 
sumptions that considerably simplify the problem. First, one works in the limit of large N^, 
in which it is readily seen [|^ that diagrams including either internal fermion-antifermion 
loops or the crossing of gluon lines at points other than their vertices are suppressed by 
combinatorial powers of Nc compared to those that do not. These simple topological con- 
sequences of the theory lead directly to the predictive power of large Nc- Second, in 1+1 
dimensions one may use the gauge freedom of QCD to choose a linear gauge in which some 
chosen component of the gluon field vanishes, and only the sole orthogonal component sur- 
vives. Then, since the gluon self-coupling term in the field strength appears as a commutator 
of field components, this term vanishes in the gauge we have selected. Consequently, gluon 
self-coupling vanishes in this gauge, and so in combination with large Nc, gluon lines are 
not permitted to cross each other, even at vertices. Moreover, ghosts are absent in linear 
gauges. It follows that the only diagrams that must be summed are "rainbow" diagrams 
for the quark mass and wave function renormalization, and "ladder" diagrams for quark- 
antiquark interactions [§. 

In 1+1 dimensions confinement is realized trivially, since the lowest-order inter-quark 
potential obtained by taking the Fourier transform of the gluon propagator (which gives 
rise to the 1/r Coulomb interaction in four dimensions) grows linearly with the inter-quark 
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separation. Although lowest-order color confinement is an automatic consequence in two 
dimensions, it is a highly nontrivial fact that the phenomenon persists in the all-orders 
Green function solutions of the 't Hooft model. 

To be specific, the Lagrangian of QCD, as in four dimensions, is 

C = -jTiF^.F^" + E (7'^(^a^ - gA^) - m,) ^„ (2.1) 

^ a 

where is the SU{Nc) gauge field with field strength F^^, defined in the usual way, and 
ipa is a Dirac fermion of bare mass rria and fiavor a. The bare coupling g not only has 
dimensions of mass in two dimensions, but scales as l/^/iVc in the large- Ai'c limit. 

The renormalization of the fermion propagator is exceptionally simple. The only modi- 
fication is a shift of the bare fermion mass by 

ml^ml^^ml-g'Nj2n. (2.2) 



Consequently, it makes good sense to describe masses in units of gy Nc/2tt, which is finite in 
the Nc ^ oo limit. The dividing line of = 1 {mlji = 0) in these units acts as a boundary 
between heavy and light quarks, as is numerically verified in Refs. [|I2|-|T^; for example, in 



p!3| it was seen that the meson decay constant approaches the standard asymptotic behavior 



oc 1/vM for M > 5 or so. It follows that g^jNc/2'K in 1-1-1 assumes a role analogous to 
Aqcd in 3-hl. 

Quantization of the theory is most convenient in axial light-cone gauge {A^ = 0), where 
light cone coordinates are defined by 

x^^x^^ (2.3) 

and analogously for other vectors. The chief advantage of this choice is that only one com- 
ponent of the Dirac algebra (7_) survives, thus effectively eliminating the need to perform 
Dirac traces. 

Upon solving for the Green function of a fermion-antifermion pair with bare masses M 
and m in this model, one obtains the bound-state eigenvalue equation 

"^(^) = + 4>T{x) - t dy 0^ (y) Pr (2.4) 



" {y-x] 



which is known as the 't Hooft equation ||^. Here the nth eigenstate (pff"^ is the meson 
wave function, the nth eigenvalue /x^ is its squared mass, and x is the fraction of the meson 
momentum's minus component (which acts, in the light-cone quantization, as a canonical 
spatial momentum component) carried by quark M. We will always label the ground state 
(the lowest mass meson) by n = 0. The principal value prescription serves to regulate the 
integrand singularity, which originates in the infrared divergence of the gluon propagator. 
This equation has a discrete spectrum of eigenvalues that increase approximately linearly for 
large n, and the wave functions vanish at the boundaries x = and 1, with the asymptotic 
behavior (p^^i^x) x^^' as a: — > 0, where 

+ ttPm cot npM = 0, (2.5) 
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and similarly as x — 1, exchanging m for M. /^m is a monotonic function of (or M^), 
increasing from zero to one as = — >■ cx). 

Also useful in this description is the full meson-quark vertex which is given by 

= /' dyrj'^iy) Ft (2.6) 
Jo [y - zy 

for all complex values of z. Indeed, except for z G [0,1], the principal value prescription is 
unnecessary. 

The decay constant jn for meson n may be computed in this framework. It is given by 



fn = Cnl\fT^-, (2.7) 



where 







rfx0„(x). (2.8) 



Strictly speaking, the r.h.s. of ( |2.7|) is also multiplied by a factor -y/iVc, but we may absorb this 
factor into the normalization of other factors by which it is multiplied in the full amplitudes; 
what is important is that the final physical amplitude has the correct Nc dependence at 
leading order. As each new quantity is calculated in this paper, we will point out the 
leading dependence on N^.^ but as a rule we suppress the explicit factors for ease of notation. 
The 't Hooft eigenf unctions 0„, for example, are 0{N^) solutions of Eq. ( |2.4|) , and so c„ is 
also 0{N^). On the other hand, light meson decay constants have the well-known behavior 
/„ oc v^iVc, and the full result Eq. ( |2.7| ), including the y/Nc, may be verified by direct 
calculation. 

The 't Hooft model wave functions (pn and are calculated by means of a standard 
numerical method called the Multhopp technique [1^,15], in which the integral equation is 



converted to an equivalent infinite-dimensional eigenvalue system, which in turn may be 
truncated after a desired number of modes to give approximate wave function solutions. 
Since the relevant formulae for unequal-mass mesons do not appear elsewhere, we present a 
summary in Appendix A. 

Whereas the eigenf unctions 4>n"^ describe the complete set of homogeneous solutions 
for the two-point Green function, the solution for 1 — 2 meson decays (the leading decay 
channels in large Nc) requires also three-point Green functions. Remarkably, the requisite 
expressions may be written entirely in terms of triple overlap integrals of the functions 
and $ without bare quark model contact-type interactions. In physical terms, this 



means that the three vertices of the diagram for the three-point Green function are resonance 
dominated, without contact contributions. Nevertheless, for the diagrams computed, it will 
prove to be computationally convenient to describe part of the full amplitude in terms of 
these contact terms. We exhibit these explicit expressions in Sec. 0, but for the moment it 
is only important to note that such expressions indeed exist. 



III. PECULIARITIES OF 1+1 DIMENSIONS 

Despite one's hopes that exact calculations in the 't Hooft model may lend insight into 
real 3+1 strong interaction physics, we emphasize that the two-dimensional universe pos- 
sesses some unique properties that must be remembered when comparing to the universe of 
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four dimensions. Therefore, the 't Hooft model may in no way be construed as any sort of 
hmiting case of real QCD, and any direct comparisons are necessarily qualitative. In other 
words, we espouse the opinion that only certain conclusions based upon our numerical stud- 
ies of 't Hooft model solutions, not the numerical results themselves, possess any validity in 
3+1 dimensions. 

The most obvious signal that 1+1 and 3+1 physics are vastly different is that the former 
does not possess the quantity of angular momentum, except in the residual form of parity]^. 
This is clear since finite rotations do not exist when there is only one spatial direction, and 
only the improper "rotation" taking — > —x^, namely parity, remains. It follows that 
't Hooft model eigenstates 0„(x) do not possess spin, but only intrinsic parity (—1)""+^ [Q. 
All of the interesting phenomenology provided by approximate spin symmetries in our world 
{i.e., the smallness of hyperfine splittings, relations between different helicity amplitudes, 
etc.) are therefore meaningless in two dimensions. 

The lack of transverse directions has important consequences for couplings in 1+1 di- 
mensions. As mentioned in the previous section, gauge couplings have dimensions of mass, 
and so such theories are super-renormalizable. Moreover, "vector" gauge bosons exist in 
1 + 1 only through their longitudinal modes. There are also different Lorentz invariants in 
1+1, since the Levi-Civita tensor e'^'^ has only two indices. The effects of these constraints 
are implicit in all the results to follow. 

The amount of Lorentz-invariant phase space is of course expected to vary between 
different spacetime dimensions D, since the measure of the phase space integrals is the D- 
dimensional volume element. However, the difference between 1+1 and 3+1 is particularly 
dramatic. To be specific, in D spacetime dimensions, the differential width for a 1 ^ 2 
decay in terms of the solid angle of either final-state particle is given by 

where |p| is the spatial momentum of either final-state particle in the rest frame of the initial 
particle of mass M, and Ai is the invariant amplitude of the process. Note particularly the 
behavior of the phase space factor |p|^~^ as the |p| = threshold is approached: For D = 4, 
the differential width vanishes with the decreasing amount of phase space available, but for 
D = 2, the differential width actually becomes singular (barring an accidental zero in the 
amplitude Ai). It follows that two-particle decay modes near threshold are enhanced in 
1+1, in stark contrast to 3+1. 

IV. THE PARTONIC CALCULATION 

Because of the small number of integrations necessary to compute phase space in 1+1 
dimensions, it is possible to perform the partonic integrals analytically in all cases of interest 



^Also, spinors retain the property of chirality, since a 75 matrix still exists in 1+1, signaling two 
inequivalent representations of the Lorentz group. 
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to us. In the 1 — 3 parton decay, one starts with 3x2 = 6 final-state momentum compo- 
nents, of which 2 are fixed by energy-momentum conservation and 3 are fixed by the on-shell 
conditions of the final-state partons; this leaves only one nontrivial integration, which can 
be done explicitly.0 

Here we consider the case of an initial quark of mass M decaying into three distinguish- 
able equal-mass quarks of mass m < M/3. The final nontrivial integral involves a small 
number of square root factors arising from the on-shell mass-energy constraints, since both 
energies and momenta appear in both the phase space and the invariant amplitude expressed 
in a given frame. Such expressions in our case integrate to the standard three kinds of com- 
plete elliptic integrals of Legendre, usually denoted by K, E, and 11. We begin by presenting 
the functional form for the phase space with constant invariant amplitude: 



$3(M; m, m, m) 



47r3M2 



where 



and 



3m , , 

^-^e[o,i]. 



(4.1) 



(4.2) 



u 



:i-e)il + e/3y 
^(l + e)(l-6/3)= 



(4.3) 



Note that, unlike the two-body phase space given by Eq. (|3.1|) , this expression does not 
diverge for finite m. However, it does possess a singularity as m ^ (e — > 0), since then 



$3 = 

The opposite limiting case e 
peculiar: 



In (—) 



(4.4) 



1, in which the three partons are produced at rest, is equally 



3v^ 
327r2M2 



[l + 0(l-e)]. 



(4.5) 



which means that phase space does not vanish in this limit. 

We now present the expressions for the inclusive partonic decay width. For definiteness, 
we attempt to describe the couplings in terms as similar to Standard Model (SM) notation 
as possible. Our labeling of partons is exhibited in Fig. 1. The decay of the heavy quark 
1 to the lighter quark 3 is assumed to couple to a vector-like weak current with vertex 



^Strictly speaking, there is also a degree of freedom from the "solid angle" of one of the final-state 
particles. However, in 1+1, this is a discrete degree of freedom. Integration of a differential width 
over this quantity gives an additional factor of 1 + 1 = 2 for Lorentz scalars and 1 + (— 1) = for 
pseudoscalars. 



7 



factor {—192/ \/^)Vzi'~i'^ (cy — caI^), carried by a gauge boson 'W of mass My/; in the SM, 
Cy = = 1/2. The couphng at the other end of the weak current, creating quark 5 and 
antiquark 4, is assumed to be the same except for the "CKM element" Vl^. is defined, as 
in the SM, by -\/25'2/8^H^i note that Gp is dimensionless in 1+1. Finally, the abbreviations 
e and u are carried over from Eqs. ( [4. 21 ) and ( |4.3| ). 




FIG. 1. Parton decay diagram for the inclusive decay. Of interest are the parton labels, as 
used in the text. 



The weak decay amplitude and width in this case are effects of orders y/W^ and Nc, 
respectively. This counting may be established in the parton diagram by observing that the 
pair (54) in Fig. 1 can occur with each of the Nc colors, but "sewing up" the partons into 
color-singlet mesons (54), (12), and (32) costs a factor of I/a/ZVc each. Finally, each color 
may occur in the loop created by 1, 3, and 2, for one more factor of Nc- It follows that the 
weak decay widthQ calculated from the parton diagram is 0{Nc). In Sec. we show that 
the hadronic calculation of the width also produces a leading factor of A^^^. 

The width is presented in two special cases. In the first, we take Mw ^ M, the usual 
four-fermion coupling assumption. This corresponds to using only the Qf^i, term in the 
numerator of the W propagator^, —i {g^^ — q^q^/M'^r) / {q^ + is)- We then find 



TC 



V 



,1/2 



X 



E(?i)-16(e/3)'(l-e/3)"'(l 



(4.6) 



The limiting cases of this expression are given by 



AGlM 



TT 



1 - - + 0(eMne) 



(4.7) 



as e 



0, and 



^The difference from the strong width, which is 0{1/Nc), is that the qqW vertices are unsuppressed 
in large while the gg-gluon vertex is 0{l/\/TV^). 

^We use unitary gauge in order to avoid the necessity of including additional charged Goldstone 
fields. 
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16G|M 



V3iv:,\'i4 



0{{l-ef) 



(4.8) 



as e ^ 1. We see that the width is finite as e — > and vanishes as e — 1. The former hmit 
shows that the width for the partonic decay of the heavy quark is given approximately by 
M times a constant, dimcnsionlcss coefficient. 

It is easy to understand the prefactor (cy — c^)^, which means that the width in the 
Mw ^ M hmit vanishes for y ± ^4 currents. The decay vertex in this hmit is of the form 
g^^ Jixjv for some quark currents J and j. Note that the only non- vanishing components of 
the metric are and g '", so the vertex involves only J-j+ and J+j-. Now, V±A currents 
correspond to the quarks being all right (left)-handed. The currents J± and j± in this chiral 
basis are just bilinears of 7±, since '^^'^^iIjr^l = ^Jij.'ipR,L- However, in 1 + 1 dimensions 
7_'0L = 7+V'R — 0, and so all currents of one chirality vanish, leading to the vanishing of 
the decay vertex. 



If we impose c 



1/4 from the beginning of the calculation, then we find that only 



pieces obtained from contraction with the terms in the W propagator survive, giving 
rise to the width 



^31^51' (e/3)'(l- 6/3)-=^/' (l + e)-V2 



X -8(6/3)^ 



+ (1 + 



;i + 6/3)^ + (6/3) (1 - e/3f 



[1 - e/3f {1 + e)E{u) 



K{u) 



+48(e/3)^(l + eV27)n(T;,ii)i, 



(4.9) 



V = 



where 

:i+6/3)(i-6) 

;i-6/3)(l + 6) 

The asymptotic expansions in this case are 

2e2 



ttM^, V3 



for e — > 0, and 



24373M^ 



l + -{l-e)+0{{l-ef) 



(4.10) 



(4.11) 



(4.12) 



for e ^ 1. Here we see that the width vanishes as e ^ but is finite as e ^ 1. Specifically, 
in the former limit the width is approximately a dimensionless constant times M^m? / My^. 



V. THE HADRONIC CALCULATION 

Matrix elements for exclusive 1—^2 meson decays are most conveniently written in 
terms of transition form factors. We identify the B meson in 1+1 as the ground state of 
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the Mm tower of resonances to which it belongs, and subsequently label it by 0. Consider 
the "tree" (T) diagram of Fig. 2, for which B = (12), where 1 is the heavy "6" quark; the 
matrix element is parameterized by 



(m(p') \qYQ\Oip)) 



{p + p'Yf+{q^) + {p — p'Yf-{.(l^) for m even, 

^^iv{p + p'Yf+{q^) + e^^^b - p'Yf-{q^) for m odd. 



(5.1) 



where = {p — p'Y, and Q and q indicate the fields of quarks with masses M and m. 
The light quark field q here refers to the daughter of the heavy quark (3 in Fig. 2), not 
the spectator quark (2 in Fig. 2), although both are taken to have mass m. The label m 
indicates the eigenvalue index of the final-state decay product meson (23) not coupled to 
the flavor-changing current. In the remainder of this section, m exclusively means this value 
and not the value of the light quark mass. 




FIG. 2. Diagram for "tree" (T) meson exclusive decay. Numbers indicate quark labels used 
in the text (except 0, which refers to the ground-state "B" meson), while letters indicate 
the eigenvalue index of meson resonances. One can also consider contact-type diagrams, in 
which the point labeled by n is not coupled to a resonance. 



Next, we reserve the label n in the T diagram for the meson resonances or contact terms 
(13) coupled to the flavor-changing current, n carries the momentum transfer g^, which 
is the kinematic variable of interest in this system; however, it proves more convenient to 
use the equivalent Lorentz- invariant quantity u = p^/q^, which indicates the fraction of 
light-cone coordinate "spatial" component of the current g_ carried by meson 0. In the 
method of calculating the matrix element (|5.1| ) developed in |]12[, one considers not mn 
directly, but rather the crossed process n 
that case, one finds u G [0, 1]: 



Om above its threshold (g > (/xq + /^m) ); in 




\ 



g2 



+ 



g2 



(5.2) 



Here and below we use the same symbol fi for the masses of heavy-light and light-light 
mesons, since from the index one can immediately tell which one is appropriate {e.g., fio is 
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heavy-light) . Since uj is obtained by solving the quadratic equation c^uP' + (/i^ — — g^)c<j + 
/Xq = 0, it should be pointed out that the branch choice used for uj does not affect the final 
numerical results for form factors or amplitudes; the two branches simply correspond to the 
two possible directions of the mesons and m in the rest frame of n. However, the branch 
chosen above turns out to greatly facilitate the numerical computations. For some values of 
below this crossed-process threshold, lo is complex, and the following expressions for the 
form factors must be computed in a different way, as discussed below. 

With the aforementioned identifications, we may express the form factors entirely in 
terms of resonance quantities, as promised in Sec. The notation and form factor ex- 
pressions we present here appear in Ref. [|I^, while the characteristic integral expression 
contained within was first obtained in ^ . The form factors are explicitly 



(5.3) 



and 



Si 



^0 



(5.4) 



where the pole residue functions An and i3„ are given by 

c„ [1 + (-1)" 



AnW 



(5.5) 



and 



\n+m+l 



(5.6) 



and the triple overlap integral F„om is defined by 



1 



l-UJ 



dv4>f{v)4>f 



23 



UJ 



1 



UJ 



1 

UJ 



l23 
m 



UJ 



UJ 



(5.7) 



We now compute the invariant matrix element for a current coupling of the form 
{—W2/ V^)V3i'y^ (cy — CA75), exactly as for the inclusive decay. Such a calculation is possi- 
ble for an arbitrary combination of V and A currents, even though we presented the matrix 
element only for a current of the bilinear = q'j'^Q, because the two currents are related 
by 7/i75 = ^fiul'^- The invariant matrix element is simply the product of a linear combination 
of the form factors determined by the current coupling, multiplied by the propagator of the 
flavor- changing current (the 'W') and finally by a factor representing the meson formed 
from the flavor- changing current. The last step amounts, via LSZ reduction of the two-point 
Green function, to the insertion of a factor of the meson decay constant. In the T diagram, 
we assign this meson the label k and quark structure (54). Using Eq. ( |2.7| ) to write the 
decay constant fk in terms of c^, we have at last 
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[Cv + CA){-lf - (Cy - Ca) [(Cy + C^)(-l)" - (cy - Ca)] 



■^^F^omiu;), (5.8) 

where the on-shell process has = /i^. The pseudoscalar parity of the ground state B has 
been taken into account in this expression. We remind the reader that in this expression 
FnOm is given by ( |5.7| ) only for n such that /i^ > (/io+/^m)^; other methods must be employed 
for smaller /x^, as described below. 

The conversion of the decay constant to in fact gives the only surviving factor of 
y/Nc in the amplitude, which means that the weak decay width is proportional to A^^^ , in 
agreement with the partonic result of Sec. |V[ This may be seen with reference to Fig. 2 by 
the usual large Nc counting arguments: The coupling of three mesons (0, m, and n in this 
case) appears with the factor while meson n is destroyed by the W current, thus 

providing a decay constant at 0{^/I^). This part of the diagram alone, which is none other 
than the form factors f±{q'^), is thus 0{N^). The creation of the meson k from the weak 
current gives the remaining factor of \^7^. Finally, the width is given by Eq. ( |3.1| ). 

The question remains what to do with the contributions from current-coupled resonances 
below the thresholds for the T diagram. The expressions listed above are inadequate because 
they assume the reality of u in the computation of contour integrals with denominators of 
the form {ut + c ± ie) , where the e arises from the Feynman prescription in the fermion 
propagators, and c represents other purely real numbers arising from the loop calculation. 
Such integrals are naturally trivial when u is real and simply lead to step functions. However, 
when uj is complex the results are rather more cumbersome (although still tractable in 
principle). One may resort instead to other methods |jl2| in order to obtain amplitudes from 
the below-threshold resonances. The approach relies on sum rules that are satisfied by 
the amplitudes, and are described in Appendix B. The upshot is that the sum rules may 
be used to describe the below-threshold amplitudes in terms of combinations of values from 
the above-threshold amplitudes, a process to which we refer as "backsolving" . 

However, backsolving has drawbacks from the practical point of view. As the number 
of resonances below threshold increases, the number of the above-threshold pole residues 
and corresponding accuracies with which these are computed must increase dramatically to 
maintain the accuracy of the below-threshold residues thus calculated. For transitions of 
the B meson to light vr's, it is known I|l2|,|l4l that the below-threshold pole residues are very 
large compared to their above-threshold fellows, and tend to alternate in sign. Clearly, a 
small uncertainty in the above-threshold calculation magnifies to a large uncertainty in the 
below-threshold residues, and the alternating sign suggests delicate cancellations among the 
computed residues, which makes the situation even worse. 

There is a much more efficient method of calculation if we are willing to abandon the 
requirement that all vertices in the calculation of Fig. 2 are resonant couplings, and allow for 
quark model-type contact terms. The calculation is based upon the observation that u as 
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defined in Eq. ( p.2| ) is real not only for decays in the crossed kinematic region > (/io+/im)^, 
but also decays in the direct decay kinematic region < < (/Uq — /im)^, where u >1. It 
therefore makes sense to redefine uo = q^/p- rather than P-/q-, so that u G [0, 1] in this 
range. One then finds 



\ 



r + K 



(5.9) 



It is convenient to define the triple overlap integral 



.i— CUJO \UJ / \l—UJ 

as well as the contact terms 



Lu Jul \UJ 



l32 



l-UJ 



(5.10) 



C2 
C3 



UJ I dv(j)Q {v 



i32 



V — UJ 



1 — UJ J v(v — Uj) 



1 — UJ Jo 



dv(j)l'{v)^ 



32 lV~UJ 
""'l-UJ 



(5.11) 



Note that the triple overlap is somewhat different from that in Eq. ( |5.7| ), both in the ar- 
guments of each wave function and the definitions (( ^.2|) and ( |5.9|) ) of uj for each case. 
Furthermore, there is some flexibility in how one expresses results containing these contact 
terms, since one can use the completeness of the 't Hooft model eigenfunctions on x G [0, 1] 
to show J2n(^n(pn{x) = 1, and from this prove identities such as 



[UJ 



Co - c,. 



(5.12) 



After a lengthy but straightforward calculation, one finds 



An{q^ 



„ 1 ~ q^ ui 



{q^Ci - [1 + (-l)'"mim3] C2 + C3} , (5.13) 



and 



E -An{(f ) f 2 _ 
. 1 - aVul ^ 



2 /,,2 

n n 



+ (l-r) (g'Ci-C; 



l + (-l 



, m+l 



Co 



(5.14) 



where 



r = 



q^/uj 



2 ' 



(5.15) 
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and the pole residue functions are given by (compare ( |5.3| ) and ( |5.4|) ) 



and 



'Sn(g^) = C„ 



l + i-l 



in+m+1 



(5.16) 



(5.17) 



Finally, the matrix element for the decay mk, which unlike Eq. ( ^.8|) holds for all such 
decays allowed by kinematics, is given by 



Mt = Ck\ - 



2 GpM. 



w 



q^- f^l 



^J^Jicv + CA){-lf - {cv - Ca) 



w 



E 



n q^ f'n 



(cv + CA)(-l)"/i„ - (Cy - CA)g -Fonm(t^) 



+ (cy - CA)g^Ci + (Cy + CA)mim3C2 



(5.18) 



where as before, = 



VI. RESULTS AND DISCUSSION 



We computed the weak decay width of a free heavy quark with masses M 
in units of (y'A/A''c/27r using Eq. 



2.28 15.00 

), the Mw — > cxD case. Likewise, we computed the hadronic 
width using the same range of heavy quark mass and a fixed light quark mass m = 0.56. The 
expressions used were Eqs. ( |5.18|) and (|3.1| ), with definitions ( ^TSf ), ( |5.9|) , ( |5.10|) , (|5.11|) , and 
with sums over all channels mk satisfying the on-shell condition Hm + f^k < A^o- Both 
widths are taken to have the same overall multiplicative factor 2G'^\VsiV^^\'^ {cy — c\Y/n. 

It is equally possible, in principle, to use ( pTB| ) instead of ( |5.18| ) and backsolve for pole 
residues An and Bn defined in ([5.5|)-(p^, or equivalently the overlap integrals -F^om, using 



the expressions in Appendix B whenever /i„ < /xq + f^m, and obtain the hadronic width in 
this way. However, as discussed in Sec. |V], this approach rapidly leads to uncontrollably 
large numerical uncertainties. Nevertheless, we were able to show in some simple cases with 
only a few backsolved residues that both methods produce the same numerical result within 
a few percent. 

It is no more difficult to consider cases other than Mw oo. For example, if one imposes 
the condition V — A condition cy = ca = 1/2, then Eq. ( ^.18|) is just as valid, but now one 
uses the partonic width ([4.91). 

Of course, the partonic width is just a single easily evaluated function of the quark 
masses. The hadronic width, on the other hand, requires first the solution of the 't Hooft 



14 



equation, which is accomphshed by means of the Multhopp technique described in Appendix 
A, repeated for as many resonances as desired. Next, the matrix elements are obtained by 
taking sums of overlap integrals over these wave functions, as in Eqs. ( |5.1CI| ) and ( |5.11| ). We 
compute the first 500 eigenvectors but include only the first 50 in our sums over resonances. 
The results change very little when more resonances are included. Finally, the amplitude 
for a given exclusive process is squared and multiplied by phase space to give the hadronic 
width. 

Clearly, such a procedure uses a significant amount of computing time, and therefore it 
is not practical to compute the hadronic width at points exceptionally finely spaced in M. 
In practice, we computed each above-threshold amplitude at values of M = 2.28 and each 
integer mass from M = 3.00 to 15.00. The significance of the lower bound is that, with the 
given light quark mass m = 0.56, this "6 quark" mass gives a ground-state "S" meson just 
above the threshold for producing two ground-state "vr" mesons, i.e., the smallest value of 
heavy quark mass unstable under hadronic weak decay. 

We then make the empirical observation that the amplitudes for exclusive processes 
AiriO —>■ mk) are smooth functions of M. We thus obtain the value of the amplitude 
at all intermediate points M by fitting to a fixed power law behavior over each interval, 
either by interpolating for values between adjacent pairs of points where the amplitudes 
were computed directly, or by extrapolating from the nearest two points if we are probing 
values of M above the process threshold but below the first explicitly computed point. In 
fact, we find that the exclusive amplitudes do not vanish at threshold and are usually]^ 
monotonically increasing functions of M (for example, see Fig. 3), although the rate of this 
increase is dependent upon the particular exclusive mode under consideration. The phase 
space is a known function of the computed meson masses, and thus the width can be reliably 
computed at any value of M in the desired range. 



^In the few exceptions to this rule, the amplitude dips slightly for values of M just above threshold, 
but thenceforth assumes monotonically increasing behavior. 
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FIG. 3. Weak decay amplitude M.t for the exclusive decay to the lowest mode — >• (m = 
0), (k = 0), as a function of heavy quark mass M, with light quark mass m = 0.56. The 

overall factor 2^2/7i CpV^iV^^^Cy — c\) in the amplitude is suppressed for convenience. 

Since the phase space in 1+1 is singular at threshold (Eq. (p.l|)), one would expect a plot 
of width r vs. M to be very ill-behaved, with dramatic singularities increasing in density 
as M increases. One would expect it to be essential to use some sort of smearing in M to 
properly test duality between this hadronic description of the F and the smooth partonic 
result. In fact, this does not appear to be the case. We refer to Fig. 4, which is our central 
result. It is obtained by interpolating each exclusive decay amplitude, as described above, 
at intervals of AM = 0.01. The remarkable result is that, after passing the first couple of 
thresholds, F appears to be a nearly smooth function in M, barely sensitive to the phase 
space singularities as each new threshold is passed. This result suggests that the effect of 
individual higher resonances is quite minimal, as one might expect in 3+1 dimensions. In 
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the 1+1 case, however, the result is all the more surprising, since now phase space near 
threshold provides a large enhancement rather than a suppression. 




FIG. 4. The full decay width for the sum of exclusive modes in the decay mk as 
a function of heavy quark mass M, with light quark mass m = 0.56. The overall factor 
SGplV^iV^^l"^ {cy — c\)^/7i in the width is suppressed for convenience. The dashed line is the 
tree-level parton result, Eq. (\i.6i). 



It is interesting to watch the width develop as more and more resonances are added. 
In Figs. 5a-5d we include exclusive channels with the lowest 1, 3, 5, and 11 thresholds, 
respectively. We now see explicitly that the full width over the range in M we consider is 
essentially produced by the first 11 channels, indicating the decreasing influence of individual 
higher resonances. The small wave in F above all the included thresholds is an artifact due 
to the interpolation routine between values of M at which the amplitudes are explicitly 
computed; its small size indicates the smoothness of the amplitudes in M and the reliability 
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of the interpolation. 
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FIG. 5a. The full decay width as a function of heavy quark mass M, with light quark mass 
m — 0.56, including only the exclusive mode with the lowest threshold value (corresponding 
to m = k = Oj. The scale is the same as in Fig. 4. 
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FIG. 56. Same as Fig. 5a, except now including the exclusive modes corresponding to the 
three lowest threshold values. 
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FIG. 5c. Same as Fig. 5a, except now including the exclusive modes corresponding to the 
five lowest threshold values. 
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FIG. 5d. Same as Fig. 5a, except now including the exclusive modes corresponding to the 
eleven lowest threshold values. Observe that this figure is almost indistinguishable from the 
full result, Fig. 4. 

Another remarkable feature of Fig. 4 is the near-perfect linearity of F for values M > 7.0. 
Suppressing the proportionality constant between F and M, Fig. 4 appears to obey the 
asymptotic form F ^ 0.514M — 0.141. This is surprisingly close to what is predicted asymp- 
totically for the partonic rate: Suppressing the same proportionality constant in Eq. ([4. 7]), 
one predicts Fp^rt = |M(1 + 0{1/M^)). 

One may ask whether the strength of the peaks in Fig. 4 is large enough that the 
mass-smeared partonic and hadronic widths nevertheless disagree. That is, local duality 
appears remarkably well satisfied, but perhaps global duality actually fails by concealing a 
large portion of F(M) in the very narrow threshold peaks. In this scenario, the apparent 
asymptotic smoothness of Fig. 4 fools us, for the density of threshold singularities increases 
with M so rapidly as to push the curve of hadronic F(M) out of agreement with Fpart(^) 
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for sufficiently large M. We now argue that this possibility does not appear to be realized, 
at least numerically. Let us smear in M over a region of size A, 1 <^ A <C M .0 Owing to the 
approximate linearity of squared meson masses in the excitation number, there are ~ M^A 
thresholds in this region, and the contribution to the smeared rate from their phase space 
near the threshold scales as M~^/^ for each.Q We observe empirically that the magnitudes 
of amplitudes ffist appearing at a given threshold mass Mthr tend to evolve approximately 
no faster than as M^^°'^. It follows that the contribution to the smeared r(M) from the 
region of width A scales approximately as M"*^ ''. For the border region, where one of the 
mesons m, k is highly excited and the other is near the ground state, the phase space is 
seen to scale as M~^, but the number of such states is only ~ MA, so again the area under 
these peaks contributes little to T{M). Finally, the phase space far above a given threshold 
scales as M~^, which means that, given the density of states for various eigenvalue indices, 
amplitudes cannot on the average grow with M above their respective thresholds faster 
than M^/^ for ^ 1 and m ^ 1, M^/^ for one of m, = 0(1) and the other ^ 1, or 
for m and k = 0(1), or else the linear behavior of r(M) will be violated. In fact, the 
amplitudes we have computed all obey these constraints. We see that the linear behavior 
observed requires a delicate balance of numbers and mass dependences of amplitudes versus 
excitation numbers, and we hope to obtain analytic arguments for this remarkable behavior 
in the future. 

What are we to conclude from this result? The 't Hooft model is exactly soluble, so 
it must be the case that the fully-dressed parton diagrams give results agreeing with the 
hadronic calculation; indeed, this is how the hadronic problem was solved in the first place. 
The partonic width computed in ( ^.6| ) represents only the Born term in an expansion in 
strong coupling g, so the addition of gluon loops is apparently necessary to bring the two 
results into agreement. The small discrepancy between the curves may have this origin, 
or it may simply be a limitation of the numerical accuracy of the calculation. However, 
it is interesting to note that the two curves appear to differ asymptotically by a constant, 
which for plots linear in M is a 1/M correction. Therefore, we suggest that this effect 
is genuine and not a numerical artifact. In Fig. 6 we superimpose on the hadronic width 
of Fig. 4 the curve Tpa.^t{M) ■ (1 + 0.15/M), and see that the fit is outstanding. From 
this result, we learn that local duality for this system is violated badly only for the first 
few resonances and very close to thresholds of higher resonances, and that 1/M effects 
appear to be only at the few percent level for M > 7. It would be very interesting to 
see explicitly what happens to the partonic inclusive width at the one- or two-loop level. 



^In practice, we use a normalized Gaussian window function with a mean of M and a variation 
of A, although the result should be independent of the particular form used. 

^This is verified from ( |3.lD and the observation that, for M ^ 1, /ig oc M. 
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FIG. 6. The full decay width of Fig. 4 compared to the tree-level parton result of Eq. ( ^.(^ 
corrected by al/M effect: Tp^^iM) ■ (1 + 0.15/M). 

One natural idea of how to improve the Born result is to replace the bare quark masses 
with the renormalized values. This would not sum all gluon corrections, but it would include 



an important subclass of them. Unfortunately, due to the result (|2.2| ), masses below 1.0 (such 
as that of our light antiquark) have imaginary renormalized values, and then our whole 
interpretation of phase space, essential for the calculation of the width, becomes ambiguous. 



VII. CONCLUSIONS 

We have calculated the nonleptonic decay width of a heavy-light meson in the context 
of the 't Hooft model as a function of the bare heavy quark mass, both for the Born term of 
the free partonic decay (which we called the "partonic width") and the full sum of allowed 
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hadronic decays (the "hadronic width"). We found that these two quantities approximately 
agree at leading order in M, with the hadronic width being slightly larger. Both quantities 
are observed to grow linearly and smoothly for large M, despite the effects of numerous phase 
space threshold singularities in the hadronic case. The slight discrepancy between hadronic 
and partonic widths is well-fit by a 1/M correction, rhadr(^) ~ rpart(Af) ■ (1 + 0.15/M). 

Assuming that the small discrepancy between the partonic and hadronic results is gen- 
uine (rather than a numerical artifact) leads one to conclude that nonleptonic heavy-light 
meson decays in 1+1 dimensions cannot be described in terms of an OPE that lacks 1/M 
corrections, and it naturally leads one to believe that the same conclusion is true in 3+1. 
Since the lowest order of the OPE is simply the naive free quark picture, this result also has 
obvious implications for the application of quark models in such decays. Another incisive 
test of quark-hadron duality in 1+1 is whether annihilation diagrams, in which the valence 
quarks in the decaying meson annihilate through a weak current, are suppressed compared 
to spectator tree diagrams (Fig. 2); these studies are well underway |]rT|, and results will be 
forthcoming shortly. 

A number of unanswered questions not addressed by this work include the effects of 
loop corrections to the Born amplitude free quark decay, the dependence of decay widths on 
the light quark mass, the effects of including finite meson strong decay widths (which are 
0(1/A^c)); the effects of identical final state quarks or mesons, multiparticle final states (also 
suppressed by powers of N^), and so on. While "two-dimensional phenomenology" cannot 
be used as a quantitative substitute for the standard four-dimensional variety, it clearly 
indicates the limitations of the standard lore. 

Note Added. An interesting recent work by Blok [^] suggests that global quark-hadron 
duality at high energies in the 't Hooft model with massless quarks may be achieved by 
including smearing through the 1 /Ai'c-suppressed widths of resonances. Our calculation, on 
the other hand, does not include finite-width effects but nevertheless achieves an effectively 
smeared result, even at relatively low M, which supports the claim of duality at leading 
order. 
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APPENDIX A: THE MULTHOPP TECHNIQUE 



This technique P^ , P^ is used to solve numerically certain singular integral equations in 
a systematic expansion of basis functions. Specifically, it is used to solve equations of the 
form 

^(x)= ['dyK{x,y)FT-^.^, (Al) 
Ja [X — y) ay 

where V'(^) = '^(^) =0- The 't Hooft equation is seen to be of this form after an integration 
by parts. Generally speaking, one maps the interval y G [a,b] to 9 E [0,7r] by a function 
linear in cos 6', and then expands in a Fourier series in 6. Equivalently, such functions 
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written directly in terms of y for each mode turn out to be the product of a common factor 
ly{l — y) times Chebyshev polynomials of the second kind, Un{y). 
In particular, we use a slight variant of the Multhopp technique for the special case that 
K{x, y) is independent of y. The variable transformation is 



1 + cos 6* 1 + cos 9' 

X = , y = , 

2 2 

in terms of which the 't Hooft equation may be written 
2 sin 9 Jo 



Expanding 



(A2) 



LMmlnt\ 



d9' 



Pr 



1 



cos 6* — cos 6*') 
(A3) 



^Mm(^) = g a^sinm^. 



(A4) 



m=l 



and using the integral (m = 0, 1 

d9' Pr - 



cos 9' — cos I 



■ cos m9' = vr 



sinm^ 
sin^ ' 



(A5) 



we are led to a series equation for the eigenvector coefficients am- 

We truncate the series after mode m = N, and evaluate both sides at the equally-spaced 
values of 9 {Multhopp angles), 



9k 



m = l...N, 



N+1' 

which are a convenient choice because of the inversion identity 

N ^ 

sin k9m sin k9n = -{N + l)Smn, 
k=i ^ 

to obtain at last the finite eigenvector system (compare Ref. [jl2|, Eqs. (A.6)-(A.7)) 



(A6) 



(A7) 



N N 



fi'^a„ 



sin 9kn sin 9km \ {M^ + m^) - (M^ - m^) cos 9^ 



sin 9k 



sin 9k 



+ mn > am- (A8) 



The series expansion ( [A4D transforms the normalization condition Jq dx (f){x)'^ = 1 into 



N N 

l = -Y,mamT. [l + (-1)'""'"] [(l -{m- nf) (l - (m + nf) 



(A9) 



m=l 



n=l 



while the phase of the eigenvectors may be chosen by noting the asymptotic forms 



N 



(f){x) -2v^ J2 i-^rmam + 0{x) 



(AlO) 



m=l 
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as X — s> 0, and 

N 

^+2v^r^^ ma^ + 0(1- x), (All) 

m=l 

as X — >■ 1. Note that the Muhhopp solution requires the eigenfunctions to vanish as square 
roots at the endpoints, in contrast to the dynamically-generated exponents of the exact 
solution given by Eq. (|2.5| ). Great care must be exercised when extracting information near 
the endpoints of the numerically-calculated wave functions. 

In terms of the eigenvector components a.^, the vertex function $ may be written for 
real x, 

(l>{x<0)= ;^(-l)"^ma^ (yr^- y/^j , (A12) 

/x{x - ^ ^ 



and 



m=l 



+7r ^ 2m 

$(X > 1) = -^^=^= ^(^rn (V+X - Vx - l) . (A13) 

Jx{x — 1) m=l 

These regions are the only ones for which we need explicit expressions; whenever x G (0, 1), 
the 't Hooft equation ( |2.4D and the definition ( ^.61 ) may be employed to rewrite $(x) in 
terms of 0(x). 

Because of the singularities in $ at x = and 1, when numerically computing overlap 
integrals such as in Eqs. (^.Tp , ( ^.101) , or ( ^.11|) , one must sample the region near the sin- 



gularities more heavily. Moreover, in contrast to the numerical Multhopp solution for $(x) 
given above, the singularities in the exact solution behave as <J)*^™'(x 0^) ~ {—x)'^'^'~^ 
and <|)*^'"(x 1+) ~ (1 — x)^™~^. Such a difference can also lead to substantial errors when 
computing overlap integrals; the answer is to increase the number of eigenvector modes in 
the wave function solution, so as to decrease the size of the region where the asymptotic 
behaviors differ. 

Finally, the quantity c^, which is the meson decay constant up to a normalization factor 
(Eq. ( |2.7| )), is simply given by 



TT 







dx(j)n{x) = — Oi. (A14) 



APPENDIX B: AMPLITUDE SUM RULES 

The asymptotic forms of 't Hooft model solutions may be used to obtain constraints on 
the form factor expressions, such as Eqs. ( |5.3| )- (|57^ ), that are derived from them, as is shown 
in 1^. There it is seen that, as ^ oo, the form factors vanish at least as fast as 
where jSm is defined in Eq. ( p.5|) . By symmetry there is also a term that falls off as 
but since M > m, we have i3m > Pm, and the slower fall-off dominates. 

It is then apparent from Eqs. (p.3|)-(p^ that the pole residue functions Aniq"^) and 
En{q^)/q^ must also vanish at least as fast as But (|5.3| ) and ( |5.4|) have a very 
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suggestive form: They explicitly display the pole-dominance nature of the large-A'^c limit, 
since Aniq"^) and Bn{q^)/q^ are expected to have no non-analytic behavior for finite g^, and 
moreover, these expressions are written very conveniently for an application of Cauchy's 
theorem. If we choose our contour such that it encloses all of the complex plane except 
for the part of the real axis where the resonance poles lie, then using Cauchy's theorem with 
the vanishing of the residue functions for \q^\ —>■ oo gives 

^l-qVfil ^l-qV^^r ^ ' 

where A„(g^) = Aniq"^) or Bn{q'^)/q^, and the sums over n are restricted to the appro- 
priate parities. Each residue function thus needs only to be computed at the value of its 
corresponding mass eigenvalue. 

This sum rule also gives us enough information, in principle, to obtain residue functions 
that are not conveniently computable in a direct fashion . Let t denote some threshold in 



n, below which it is unwieldy to compute the residue functions directly. Then, if we choose 



at least t values of qf, i = 1 . . .t above threshold, we obtain a solvable t x t linear system of 
equations: 



-^l-oVu^ ^l-o2/u2 ^l-oVu^' ^ ^ 

n<t y« / r^n n>t / h^n n Hi I H'n 

Adding additional qf simply overconstrains the system and provides consistency checks. We 
refer to this approach in the text as "backsolving" , to indicate that it is an indirect method 
of solution. 

We also note one other sum rule [p|,p!2|, which is simply derived by taking the limit 
— >■ oo in ( |BT| ) and recalling that the form factor expressions vanish in this limit: 

^/iX(/in) = 0. (B3) 
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Figure Captions 



FIG. 1. Parton decay diagram for the inclusive decay. Of interest are the parton labels, 
as used in the text. 

FIG. 2. Diagram for "tree" (T) meson exclusive decay. Numbers indicate quark labels 
used in the text (except 0, which refers to the ground-state "B" meson), while letters indicate 
the eigenvalue index of meson resonances. One can also consider contact-type diagrams, in 
which the point labeled by n is not coupled to a resonance. 

FIG. 3. Weak decay amplitude Mt for the exclusive decay to the lowest mode — *• 
(m = 0), (k = 0), as a function of heavy quark mass M, with light quark mass m = 0.56. 
The overall factor 2^2/71 CpV^iV^^^Cy — c\) in the amplitude is suppressed for convenience. 

FIG. 4. The full decay width for the sum of exclusive modes in the decay mk as 
a function of heavy quark mass M, with light quark mass m = 0.56. The overall factor 
8Gp| V3iV"4*5p(cy — c^)^/vr in the width is suppressed for convenience. The dashed line is the 
tree-level parton result, Eq. ([4.61). 

FIG. 5a. The full decay width as a function of heavy quark mass M, with light quark mass 
m = 0.56, including only the exclusive mode with the lowest threshold value (corresponding 
to m = k = 0). The scale is the same as in Fig. 4. 

FIG. 56. Same as Fig. 5a, except now including the exclusive modes corresponding to 
the three lowest threshold values. 

FIG. 5c. Same as Fig. 5a, except now including the exclusive modes corresponding to 
the five lowest threshold values. 

FIG. 5d. Same as Fig. 5a, except now including the exclusive modes corresponding to 
the eleven lowest threshold values. Observe that this figure is almost indistinguishable from 
the full result. Fig. 4. 

FIG. 6. The full decay width of Fig. 4 compared to the tree-level parton result of Eq. ( |4.6|) 
corrected by a 1/M effect: Tp^,tiM) ■ (1 + 0.15/M). 
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